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Abstract. We study a 2-scale version of the Landau-Lifshitz system of ferro- 
magnetism, introduced by Starynkevitch to modelize hysteresis: the response of 
the magnetization is fast compared to a slowly varying applied magnetic field. 
Taking the exchange term into account, in space dimension 3, we prove that, un- 
der some natural stability assumption on the equilibria of the system, the strong 
solutions follow the dynamics of these equilibria. We also give explicit exam- 
ples of relevant equilibria and exterior magnetic fields, when the ferromagnetic 
medium occupies some ellipsoidal domain. 
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1 Introduction 

Hysteresis is a widely studied, yet not completely understood phenomenon. It 
has played a role from the very beginning of the works on magnetism. Lord 
Rayleigh [9] proposed a model for ferromagnetic hysteresis in 1887, while the 
most achieved micromagnetism theory goes back to Landau and Lifshitz, in 
1935 (see [ ]). 

In [12], Visintin gives many historical references, underlines the links between 
several forms of hysteresis (in particular, from plasticity, and from ferromag- 
netism), and how it is related to phase transitions. He performs a mathematical 
study of the so-called hysteresis operators, including the most famous one, due 
to Preisach. 

Recently in [3], Carbou, Effendiev and Fabrie have proved the existence of 
strong solutions to a model of ferromagnetic hysteresis due to Effendiev. 

In this paper, we rather investigate properties of a two-scale model intro- 
duced by Starynkevitch in [ ]. This model describes the dynamics obtained 
when some exterior magnetic field is applied to the ferromagnetic material under 
consideration, while the response of the magnetization occurs on a much shorter 
time scale (say, denoted by e > 0). Mathematically, such models, associated to 
ordinary differential equations, had been studied in the nonstandard analysis 
framework, leading to "canard cycles" (see [6]). Considering a Landau-Lifshitz 
model in space dimension (thus, an ODE), Starynkevitch studies the possible 
equilibria of the system, and the asymptotic behavior of the solutions (as the 
above mentioned parameter e goes to zero) when the exterior magnetic field 
slowly varies. 

Our aim is to extend Starynkevitch's approach to the Landau-Lifshitz model 
in space dimension three, taking exchange term into account. This means, giving 
the asymptotic description of solutions to the slow-fast corresponding system 
of partial differential equations. Here, we prove such a result away from the 
bifurcation points of hysteresis loops. More precisely, assuming that the system 
(described by its magnetization) possesses at each time t some stable equilibrium 
w eq (i), and is submitted to some slowly varying exterior magnetic field, we 
show that the magnetization follows the dynamics of TO oq . We also give explicit 
examples (for ellipsoidal domains) of relevant equilibria and exterior magnetic 
fields. 



2 



2 Statement of the results 



The initial and boundary value problem associated to the 2 scale Landau-Lifshitz 
equation considered reads: 

!ed t m e = m £ A h £ T — am £ A (m e A h £ T ), for t ^ 0, x e fi, 

The unknown is the magnetization m e , function of the time variable t J? and 
of the space variable x € CI, with values in the sphere S 2 CM 3 . The domain 
occupied by the ferromagnet is a subset of R 3 . Furthermore, h e T = hx(t, m £ (t)), 
where the total magnetic held hr is defined by 

(2.2) h T (t,m) = Am + h d (m) + h eKt (t). 

Here, the first term Am is the "exchange term" , which tends to impose a con- 
stant magnetization (domains where magnetization is constant are called "Weiss 
domains"), and Am denotes the extension of Am by out of f2. The second 
term, yielding spatial variations of the magnetization, is the "demagnetizing 
field" /id(m), which results from a quasi-stationary approximation of Maxwell's 
equations; it is defined (at least, for m € L 2 (fl, M 3 ), as an element of L 2 (R 3 , R 3 )) 
by 

curl /id (m) = and div (hd(m) + m) = in R 3 . 

Classical properties of the mapping m i— > hd(m) are recalled in Section 3.1. The 
third term, h ox t, denotes some given exterior field, which is assumed to depend 
on time (and possibly on space). The positive constant a is some damping 
coefficient, which appears in the model when passing from a microscopic to a 
macroscopic description. The small parameter e > expresses the fact that, 
while the exterior field h ex t depends on i, and has time variations at scale 1, 
the magnetization m E essentially depends on i/e, and thus has variations at the 
much more rapid scale e. 

Throughout this paper, for any s 6 N, we denote by H s (fl) the usual Sobolev 
space of functions with values in some vector space R , whereas H s (£l, S 2 ) is 
the Sobolev space of functions with values in the sphere S 2 (which is not a 
vector space), 

H s (fl, S 2 ) = {m e H S (Q) | \m\ = 1 almost everywhere}. 

Finally, for s ^ 2, Hfj(n,) denotes the subspace of functions in H s (il) with 
homogeneous von Neumann boundary condition, 

HUn) = {meH s (n) \d„m lan =0}, 

and 

H s N (n, s 2 ) = H s {n, s 2 ) n h s n (q.). 

All these spaces (even if not vector spaces) inherit the (metric) topology given 
by the usual norm on H s (il). 
We prove the following 
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Theorem 2.1. Let ft be an open and bounded subset ofM 3 , with smooth bound- 
ary. Let T > 0, and h cxt G C 1 ([0, T], C°°(R 3 )), bounded with bounded deriva- 
tives. Assume that there existm eq G C 1 ([0, T], i?jy(Q, S 2 )) andmo G Hff(n,S 2 ) 
such that 

(i) for all t G [0,T], m oq (io) *s an equilibrium for 



(2.3) 9tm = m A hr(to, m) — am A (m A hx{to, m) 

(see (3.1)); 

(ii) the solution hq to the initial and boundary value problem 




is global (n G C([0, oo), H^n, S 2 ))), with VAn G i 2 ((0, oo) x O), and n (t) 
converges in H 2 (fl), as t goes to oo, towards m cq (0); 

(Hi) the linearized operator £(m eq ) given by (4.31) has the following dissipation 
property: 

(2-5) 

there exist Cy m > and n > such that, 

for all 5 G C([0,T],H°°(Q)) with |m cq +S\ = 1 and d v 8\ Ba = d v M\ gn = 0, 

sup ||^(t)||fl" 2 (n) ^ V implies: 
te[a,T] 

Vt G [0, T\, (£(t,m oq (t))<5(t) | *(*)) < -Ciinll^Wll^cn). 

Then, there is e$ > suc/i i/iat, /or a/Ze G (0, Eq), the solution m e to (2.1) exists 
up to timeT (m e G C([0, T], H^(fl, S 2 ))), and converges in L 2 ((0, T), H 2 (VL)) n 
C([t, T], H 2 (Q)) towards m cq as s goes to zero, for all t G (0,T). 

To prove Theorem 2.1, we first show that m £ converges to m oq (0) within 
an initial layer of size t e = Celn(l/e). This is achieved via classical energy 
estimates (in H 2 ), carefully controlling the dependence upon e -more technically 
speaking, the quasilinear and elliptic degenerate system of PDE's in (2.1) is first 
converted into a perturbation of some linear, strongly elliptic system, yielding 
the usual smooting properties, and a Galerkine approximation is used. In a 
second step, we prove that m e converges towards m oq on the whole time interval 
[t e , T\. This amount to proving of long-time existence and return to equilibrium 
result for small initial data. Toward this end, we use again energy estimates, 
together with the stability assumption (2.5). 

Figure 1 illustrates this corresponding asymptotic behaviour. 

The above assumptions on the equilibrium m eq are discussed in Section 3.2 
below. In particular, Assumption (ii) in Theorem 2.1 may be understood as a 
choice of 'prepared' data tuq allowing to deal with the initial layer (0, celn(l/e)). 
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Figure 1: Dynamics of the magnetization away from bifurcation points. 

Figure 1: Dynamics of the magnetization away from bifurcation points. 

The dissipation property (2.5) expresses, for all to £ [0, T], the stability of the 
linearization around m cq (io) of (2.1), with e = 1 and with h CKt replaced with 
h e xt(to), independent of time. This is a strong assumption, which ensures global 
existence of the solutions to the corresponding Landau-Lifshitz equation, for 
initial data close to m cq (to)- 

Proposition 2.2. Let Q be an open and bounded subset of R 3 , with smooth 
boundary. Consider an exterior magnetic field h cxt £ C°°(R 3 ) (independent 
of time) bounded with bounded derivatives. Assume that there exists m eq € 
i5/^,(f2, 5 2 ) (independent of time) satisfying the equilibrium condition 

(2.6) m cq A (Ara oq + h d (m cq ) + h cxt ) = on ft, 

as well as the stability condition 

there exist Cn n > and rj > such that, 

for all 5 E H°°(ft) with \m cq + S\ = 1 and d v 5\ an = d v M\ an = 0, 
( 2 - 7 ) ||5(t)||ff2(n) ^ V implies: 

(c(0,m cq )6\s) H ^ -C lin \\S\\ 2 H2m , 
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for the linearized operator £(0,m eq ) given by (4.31) (with m cq (0) and /i e xt(0) 
replaced with m cq and h ex t, respectively). 

Then, there exists rjo > such that, for all mo <E Hjj(n,S 2 ) satisfying 

1 1 TOO - TO e q||jf2(r2) < 

i/ie solution n to the initial and boundary value problem 

{d t n — n A fyr(Q, n) — an A (n A /ir(0, n) j, 
n u=a = mo, 

is global (n e C([0, oo), flj^fi, S 2 )) J, wtfi VA(n - m cq ) € £ 2 ((0, oo) x O), and 
n(i) converges in H 2 (fl), as t goes to oo, towards m eq . 

In the case of m eq (0) constant over fi, Proposition 2.2 expresses that in 
Theorem 2.1, assumptions (i) and (hi) imply assumption (ii), so that we get: 

Corollary 2.3. Let ft be an open and bounded subset of R 3 , with smooth 
boundary. Let T > 0, and h ext £ C 1 ([0, T], C°°(IR 3 )), bounded with bounded 
derivatives. Assume that there exist m cq e C 1 ([0, T], Hfj(Q, S 2 )) satisfying as- 
sumptions (i) and (Hi) from Theorem 2.1. Assume furthermore that m oq (0) is 
constant over Q. 

Then, there exist r]o,£o > such that, for all mo € Hfj(n, S 2 ) such that 

\\m - m cq (0)\\ H 2 {n) ^ rjo, 

and for all e € (0,£o)> the solution m e to (2.1) exists up to time T (m e £ 
C{[0,T],H 2 N {n,S 2 ))), and converges in L 2 ((0,T), H 2 {Q))nC([t,T], H 2 (n)) to- 
wards m eq as e goes to zero, for all t € (0, T). 

In Lemma 3.5 below, we give examples (in ellipsoidal domains) of equilibria 
TO eq satisfying the assumptions of Corollary 2.3. 



3 Preliminaries 

3.1 Some functional analysis 

In this section, we recall some functional analysis results useful in the sequel. 
The first of them deals with the continuity properties of the demagnetizat- 
ing field operator hd, immediately deduced from the Fourier representation 

M^)(0 = -M(£)) j4: 

Lemma 3.1 (hd properties). Let f2 be an open subset ofM 3 . For all s inN and 

u in H S (Q), one has 

\\hd(u)\\ H s (R 3) < |M|jJ.(fi). 
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Furthermore, for all v in L 2 



we have 



{h d (u) | u)z,2 ( n) 



(u | hd(v)) L '(ci)- 



In addition to the usual Sobolev embeddings, we recall the following esti- 
mate, which results from the coercivity of the operator A = 1 — A, with domain 
D(A) = {m G H 2 (tt) | d v m\ 3U = 0} (see for example [5]) 

Lemma 3.2. Let fl be a smooth bounded open set in M 3 . There exists a constant 
C > such that for all u in H^(fl) one has 



In the sequel, we will need the following definition. 

Definition 3.3. Let il be a smooth bounded open set in K 3 . For k e N* , let P k 
be the L 2 (f2) -orthogonal projection onto Vk, the vector space spanned by the first 
k eigenfunctions of A = 1— A, with domain D(A) — {m G H 2 (fl) | d v m\ Bn = 0}. 

The family of operators (-Pfc)fceN satisfies useful properties: 

Lemma 3.4. The following properties are true. 

(i) Mk G N* ; Vu G D(A), AP k u = P k Au, 

(ii) Vfc G N* 7 Vs G N, Vu G H S {Q), P k u G H S {Q) 
(and P k u G Hfj(£l) when s ^ 2), 

(Hi) Vs G N, lim ||(1 — Pk)u\\ii s m) — f or a ^ u G H S (Q,) when s = 0, 1, and 
for all u G H^(£l) when s 2. 

Proof, (i) For all u in I? (A), fc in N*, one has 



with (V'i)iGN the L 2 -orthonormal basis of the eigenvectors of A associated to the 
eigenvalues (Xi)ieN- Then, using the vanishing Neumann boundary conditions, 



NU»(n) < C (||u||i a(n) + HAull^n)) 



P k Au = Y,( Au I *l>j)v>4>: 





so that 



k 




Point (ii) follows from the regularity properties of the family (i/>j)j e N- 

Point (iii) is a consequence of the fact that u M> ( 52jLiO- + ^'j)\( u I V'j )i 2 1 2 
provides a norm equivalent to the usual one on H s (fl). 



1/2 



□ 
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3.2 About equilibria 

Global solutions and equilibria. In [1, Th. 4.3], in the case of ellip- 
soidal domains 57 C M 3 and under a smallness assumption (on H/iextlU 00 and 
||Amo||i2), Alouges and Beauchard construct global smooth solutions to (2.1). 
Furthermore, these solutions satisfy 

VT > 0, ||Am(T)||| 2(n) + G [ ||VAro||£ a(n) < ||Am |U 2( n), 

Jo 

so that VArn belongs to T 2 ((0,oo) x 57). This is a part of our assumptions 
on the equilibrium m eq , when requiring the existence of the global solution Uq. 
Saying that ro eq (to) is an equilibrium for (2.3) means 



(3.1) 



m cq (t ) A h T (t ,m cq (to)) = 0, 
d u m cq (t ) lgn = 0, 



and requiring H 2 convergence of rio(t) towards m cq (0) as t goes to oo implies 
that m oq (0) is an equilibrium for (2.3) with to = 0. 

Energy minimization. It is worth noting that energy decay occurs along 
the evolution of no(i), so that one may hope at least H 1 convergence of no(t) 
towards some local minimum of the energy, as t goes to oo. To the Landau- 
Lifshitz system (2.1) is associated the energy 

£(t,m) = - / |Vm| 2 - - / m - h A (m) - [ m - h^t), 



2 in' ' 2_ 
and when m is solution to (2.1), we have 

^£(t,m(t)) = ~\\m(t) Ah T (t,m(t))\\l 2(n) ~ f ™{t) ■ d t h ext (t). 

Since the exterior magnetic field does not depend on time during the evolution 
of no , we get 

— E (t, n {t)) = -a\\n (t) A h T (t,n (t))\\l Hn) . 

In the case of ellipsoidal domains, special configurations are available. See 
[8], and references therein: there exists a real 3x3 definite positive diagonaliz- 
able matrix D giving the demagnetizing field resulting from any magnetization 
constant constant over 57: 

Vd £ R 3 , h d (v) ]n = -Dv. 

Hence, if u € S 2 is an eigenvector of D associated to the eigenvalue d > 0, and 
if the exterior magnetic field is h cxt = Ait for some A > (or h CKt (x) — \\{x)u 



8 



for some x € C^°(]R 3 , [0, 1]) to get a spatially localized field), then the system 
possesses two explicit equilibria rn^ q and m~ q : 

(3.2) m± = ±u. 

One easily computes the energy associated to perturbations of these equilibria: 
for all S € H%(n,R 3 ) such that \mf q + S\ = 1 a.e., 

S(mf <l + 5)-£( m t q } = l f \V5\ 2 -\I 5-h A {S) + h±\-d) ( \5\\ 



The first two terms are non- negative, so that for A large enough (A > of), rn^ q is 
a global minimum of £; but for A small, it may fail to be even a local minimum. 
Concerning m~ q , for all A > 0, if d is the largest eigenvalue of D, and 6 is 

constant in space, then the difference of energies above is less than f |<5| 2 , 

2 J n 

thus negative, whereas for S with large variations, the gradient term dominates, 
and the energy difference becomes positive. Hence, m~ q is always a saddle point 
for£. 

The dissipation property (2.5). We have the following lemma, the proof of 
which is postponed to Section 6.1: 

Lemma 3.5. For A > large enough, the equilibrium rri£ q from (3.2) satisfies 
the dissipation property (2.5) (for some constant Ci; n depending on X). 

For m~ q , it is shown in Section 6.1 that for A large, we have on the contrary: 

Lemma 3.6. For A > large enough, there exist C = C'(a, A) > and 
n = r)(a,X) > such that, for all 5 € C([0, T],H°°(Q.)) with |m cq + 8\ = 
1 and 9i,<5| on = <9„A(5| on = ; when \\5\\H^(n) ^ n, we have: 

Vt € [0,T], (r(t,m- (*))*(<) | *(*)) > C|Hl!p(Q)- 

4 Proof of Theorem 2.1 

First, consider the solution no to the Cauchy problem 

d t n = n A ft/r(0, n ) - an A (no A fyr(0, no)) , 
9 v n \ 3O = 0, 
no| 4=0 = mo, 

and define n e by 

Vt^O, rc £ (t) = n (*/£)• 
Then, n £ £ C b ([0,oo),H 2 (fl)) (with VAn £ € i 2 ((0,oo) x Q), and we know that 

(4.1) i e /e— >oo =^ n £ (t E ) — >m oq (0) in # 2 (Q). 
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Next, as in [4], we observe that for smooth functions m with constant mod- 
ulus (w.r.t. x), one has m ■ Am = — 2|Vm| 2 , so that smooth solutions to (2.1) 
cquivalcntly satisfy 

{edtrrf — aAm £ = J-(t, m £ ), 
mf t=0 =m , 

where 

(4.3) T{t, m) — m A hxit, m) + a\Vm\ 2 m — am A (m A (h^{rn) + h cx t{t))^j . 



Furthermore, smooth (L^°H 2 ) solutions to (4.2) issued from mo with constant 
modulus, equal to one, are shown to keep the same modulus for all time, (due to 
uniqueness of the solution a — \m\ 2 to ed t a = aAa + 2a| Vw| 2 (a — 1), d v a^ an = 0, 
a| t=0 = 1). We thus solve (4.2) in the Banach space C([0, T], H%(£1)), and 
deduce from this conservation that the solution actually belongs to the space 
C([0,T],i^(ft,S 2 )). 

It is worth noting that (2.1) is an initial and boundary value problem for 
some quasilinear and parabolic degenerate operator, which is seen in (4.2) as a 
perturbation of a linear and strongly parabolic one. 

Standard energy estimates ensure local-in-time existence and uniqueness of 
solutions continuous in time, with values in H 2 (Q)) (with an existence time 
depending on e): see for example [1] or [4]. By the usual continuation argument, 
we simply need to bound the H 2 norm of m 6 to ensure existence up to time T . 
Actually, we shall prove convergence (as e goes to zero) at the same time, via 
energy estimates. 

We hrst show that, after some time t e of the form t £ — Celn(l/e), m e and rf 
are close: supj 0j t e ] \\m £ — rf\\H^(Q,) goes to zero with e; thus, for e small enough, 
m e (t £ ) is as close (in H 2 (fl)) to m cq (0) as desired. We then use the stability 
property of m eq (t) to show that m e {t) stays close to it, for t E [t £ ,T]. 

4.1 First step: the initial layer [Q,t e ] 
4.1.1 Galerkine scheme 

For k £ N*, let Pk be the L 2 (f2)- orthogonal projection onto Vk, the vector space 
spanned by the first k eigenfunctions of A = 1 — A, with domain D(A) — {m £ 

l \an 

(4.2) by: 



H 2 (il) | <9„mu n = 0}, as in Definition 3.3. Define a Galerkine approximation of 



f ed t m E k - aAm| = PkF(t, m%), 
\ ml u=o =P k m . 

The projection n| = Pkrf also satisfies 

ed t nl - aA.nl =P k F(Q, n%) + a[P k ,A]n e + P k (>(0, n £ ) - JF(0, n|)) 
=P k T(0,n%) + P k [P k ,T(0,-)}(n"), 
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since for u £ D(A), PkAu = APf-u, according to Lemma 3.4. 

Now, perform energy estimates (in L 2 ) for ip e k = m| — n|, solution to 



(4.5) | 
4.1.2 L 2 estimates 



Take the scalar product (in L 2 (fl)) of <p\ with the first equation in (4.5) to get 
£ d 



2dt 



(llVfclli"(n)) + a||V^|||» ( n) = h + h + h + h, 



with 



h = (vl I ™l A h T (t,m%) - n% A /i T (0,n|) 



L 2 (Q) 



L 2 (n) 



I 2 = a(^| | |Vm||^m| - \Vn%\ J n% 

h = -a(yt | m% A (mf A (/i d (m|) + h cxt (t)) 

- n% A (n% A (h d (n%) + h cxt (0)) 

h=(if% | [P fe ,-F(O,-)]0 



L 2 (f2) 



L 2 (n) 



Estimating Ii . Decompose m| = n| + tp £ k . For all y>, h £ R 3 , ip ■ (cp A h) = 0, 
so that 



7j =(ip e k | nf A (A(n| + <p%) + h d {n% + <p%) + h cxt (t)^j 
- n£ A (An% + h d (n e k ) + h cxt (0) 



L 2 (n) 

= | n% A (A<p% + Mv>f))) + (ft I n k A (he*t(t) - h cxt (0))^ j () 

Using the continuity of h d on L 2 (S1), we get, for some constant C depending on 
||<9tfrext||L~ x and ||no||L~((o,oo),L 2 (Q)): 

(4-6) h ^CM\\ L 2 (n) (M\\H H n)+t). 

Estimating I2. Write 

|Vm|| 2 mI - |Vn|| 2 n| = (|Vm|| 2 - |Vn|| 2 )m| + \Vn%\\l 

= (V(2n| + p£) • V^)K + ¥>1) + |Vn^|V*. 
Then, use Sobolev's embeddings, such as 

(fl I (Vn| • V<p|)ri£j < ||<y5fc||L-(0)||Vn|.|| L 2 (n) ||V(^||| L 4 (n) ||n||| L 4 (0) 

^ IIVfcllfl-^^llnlllH-i^llV^llljji^llnllljj^n), 
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L 2 (0) 



L 2 (0) 



and 

(cpl | iVnUVl) < ll^lli~ ( o)l|V^|| 2 L2(0) < ||^||^ ( n)l|nl||^ 1(n) 
to get the estimate 

(4-7) I 2 «S C||(p||| ff2( o)(|k||| ff2 (n) + HvSII^(n)). 

for some constant C depending on ||no||L<=o((o l oo),_ff 1 (J2))- 
Estimating I3 . As for ii , cancellations allow to write 
h=- a(v% I n% A (n| A h d (tp%) + <p% A /i d (n|)J + n| A ^| A h cxt {t) 
~ aUp% I n| A (n% A (/i ext (t) - h ext (0)) 

Boundedness of on LP for finite p provides the bounds 

A (n| A h d {tp e k ))\\ L 2 {n) ||n||||6 (0) ||(y5||| L 6 (0) , 

and 

IK A (Vl A M n fc))ll£»(n) < l|n|||i6 ( n}||vllU«(n}- 
The above L 6 norms are controlled by 77 1 norms. Thus, for some constant C 
depending on ||/i ext ||i~, ||<9 t /i e xt|U~ and IKHl^Coo),// 1 ^)) 1 

(4-8) h^C\\yl\\ LHn) (\\ipl\\ H2m +t). 

Estimating 7 4 . Setting r| = ||[P fe ,J"(0, ')]( n f )lli,2(fi)! we have: 
(4.9) 

h < ll^llli^m + and r| — > in L^O, T) for all T > 0, with e fixed. 



This is a consequence of the following lemma, the proof of which is postponed 
to Section 6.2. 

Lemma 4.1. For all T > and n G C([0, T], ff^(fi)) n i 2 ((0, T), H 3 (fl)), 
[P fe ,J-(0,-)](n) — ► J «L 2 ((0,T),i/ 1 (r!)). 

fc— f OO 

Gathering L 2 estimates. Adding (4.6), (4.7), (4.8) and (4.9), we get 
~(ll^lll.( n ,)+ "II Vtflli.<n) 



(*+ll^ll 



'fc )) 



for some constant depending on the quantities ||ftext||i?° ) H^t^extHi? and 
IKo|U~((o,oo), m(n)), and r| from (4.9). 
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4.1.3 H 2 estimates 

Next, take the scalar product in L 2 (il) of A 2 ^! with the first equation in (4.5) 
to get 



2dt 



(||A^ fe || 2 2(f2) ) + a||VAvj|||i a(n) = Ih + Ih + Ih + Ih, 



with 



Ih = (AVl I m s k A h T (t,m%) - n% A h T {0,n%)) 
II 2 = a(&\% | \Vm%\ 2 m% - |Vn|| 2 nf 



£ 2 (f2) 



II. 



L 2 (n) 

a(AV* I m% A (m| A (h d (m%) + h ext {t)) 

- n% A (n% A (h d (n%) + h cxt (0)) 

Ih= (AVl I [i\,.F(0, •)](«") 



L 2 (f2) 



L 2 (fi) 



Estimating II\- Split 

Ih=Ih A + Ih, 2 + II h3 , 

with 

//!,! = (a 2 ^ I m| A (A^| + hd(^))) 
Ih a = (A 2 ^| I ip%Ah T (t,n%) 
7J 1)3 = (a 2 ^ I nf A (/wW - ^cxt(O))) 
The first term is written 
II 1A = (AV fe I n|A(Av»|+/»d(v*))) £ +(AVl I PfcA(Av»|+ftd(^)) yi _, i(ji 



L 2 (f2) 
L 2 (Q)' 

L 2 (S~2) 



Integrating by parts, 
(AVI |n|A(A^ + ha(^)) 



VA^| | Vnl A (A^ k + hd(<p%)) + n| A V/i d (¥>l) 



L 2 (0) 



< ^IIVA^IH^^ + ^||Vn| A (A^ + ftj^)) + n% A V/i^DH^o) 

< r/llVA^II 2 ,^ + C„ (||n|||^ 2(0) + HVAnlHlajn)) ||^|| ff 2 ( a), 
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for all T) > 0, for some (large) constant C v , using Sobolev's inequalities. ^From 
this, we deduce that for all r\ > 0, there exists C v > 0, depending only on 

lho||i~((o,oo),H 2 (n)), such that 



(4.11) /L 2 (n) 

r?||VA^|||2 2(n) + C„ (l + ||VAn£||f 2(n) ) Mf H 2 {ny 

Integrating by parts again, 
(A VI I <P% A (A^| + h d { V %)) 



L 2 (Q) 

(vA<p% | V<p% A (A<p% + h d (ip%)) + <p% A V/m(^) 



L 2 (fi) 



< llVA^III^^IIV^IIIioo^llA^III^^) 

+ l|VA^.|| L2(f2) (||Vvj|||£*(n)||/»d(Vfc))IU«(n) + ll^fclU«(n)l|V^d(^))IU*(n)) 

< ||VA^||| i2(n) (||AV^||U 2(n ) + ||V^||U 2( o)) yUnnn) 
+ ||VA^|| L2(f2) ||^||^ 2(n) , 

using ||V/id(vl)IU 4 (n) ^ ||frd(v!)llff 2 (Q) ^ ll</?|l|ff 2 (n)- Hence, there exists an 
absolute constant C > 0, and for all ij > 0, there exists C,, > such that 

(AVfcl^fA(A^f + M^t))) r2/0 ,< 

(4.12) V /L 2 (f2) 

(v + C|^||| ff 2 ( o)) ||VA<^|||2 2(n) + (7,11^111^(0). 

Summing up (4.11) and (4.12), one gets C > and, for all ry > 0, a constant 
C ri >0 (depending on |K||L~((o : oo)„F/ 2 (n))) such that 

Ih,i «S (?y+C||^|| ff 2(n))||VA^|||2 2(n) 

+ C 7) (l + ||VAn||| i2(n) + ||^|||^ (n) ) ||^|||^ 2( n). 

The second term is 

Ih.2 = - (vA<p% | Vifl A (An| + h d (n%) + W*))) 

- (VA^ | pi A (VAnl + Vh d (tp%) + v/l ext(*))) ■ 

Using Sobolev's inequalities again, we have, for all -q > 0, a constant G n (de- 
pending on ||n ||L^((o,oo),_H- 2 (o)) and ||/i xt||L t °°J such that 

-(VA^| | V<p%A(An% + h d (n%) + h ext (t))) ^ ||VA^||| L2(a) x 

x (l|VvII|z,4(n)||Anl + h A {n £ k )\\ LHn) + \\V(p%\\ L 2(n)\\heKt(t)\\ L <*>(n)) 
(4.14) <a||VA^||| i2(n) ||^||| H2(n) 

(l|VA?4|| L 2 (n) + IKIIff 2 (fi) + ll^cxt(<)|U~(0)) 



< »7||VA^||i 2(n) +C n (l+ ||VAn|||| 2(0) ) M\ 



e ||2 
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In the same way, for all 77 > 0, there is C v > (depending on ||^o||L°°((o,oo),ff 2 (si)) 
and ||V/iext||i™ x ) such that 



(4.15) 



(VApf I <p% A (VAn| + Vh d (<p%) + V/iextW) < 



L 2 (^) 



< »7l|VA^|||, (n) + C v (l + ||VAn^||| 2(n) ) ||^I||^ (n) . 



Summing up (4.14) and (4.15), we get, for all i] > 0, a constant > (de- 
pending on \\h cxt \\ L ^, ||V/i 0Xt ||L t ~ and ||nolU-((o,oo),ff 2 (Q))) such that 

(4.16) II 1>2 »/||VA^||i a(n) + C v (l + ||VAnf || 2 2(fi) ) ||^|||^ (n) . 
The third term is 

//i j8 = (a^I I ng A (^(i) - WO))) • 

Integrating two times by parts, it is easily estimated, thanks to a constant C 
depending on \\d t h cxt \\ L ^ w ^ , as 

(4-17) Ih,^Ct\\n%\\ H * { n)M\\H>m. 

This gives finally, summing up (4.13), (4.16) and (4.17): there is C > 0, 
and for all i] > 0, there is C v > (depending on 77, ||no||L°°((o,oo),_f/ 2 (n))) 
H^cxt II z.f° wi-°° and \\ 9 th e xt\\ L °° w 2,°°), such that 

(4.18) 

Ih < (v + C\\ <pi\\ H 2 (n) ) II V A I, (n) 

+ (l + ||VAn£|| 2 2(n) + |K||!p (n) ) M\\ 2 m{n) +*lkIIU=(n))- 

Estimating II2. Split 
with 

7/2,! = (A VI I (|Vm|| 2 - |Vn|| 2 )(n| + p£) 
= (A VI I V^| • (V^| + 2Vn|)m 



L 2 (Q) 



L 2 (f2) 



I/ 2)2 = (AVI I Wnl\ 2 (mt-nl)) L ^ = (aV| I V<V£ /; 
Then, using in particular the Sobolev inequality from Lemma 3.2 
||Vn|||i- ( n) < l|Vn||| i2(0) + ||VAn||| L3( n), 
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we get: 

Ihs = - (VA^| | A^ k ■ (V^£ + 2Vn|)(n| + 

- (VA^| | • (Apg + 2An|)(n| + ^)) i2(n) 

- (VA^| | • (V^| + 2Vn|)(Vn£ + V^f))^ 
sc ||VA V ||| L 2 (n) ||A^||| i2(n) x 

x (l|V^|||ioo ( n) + 2||Vn|||£oo ( n)) (KHl~(^) + Mh^(n)) 
+ l|VA ¥ p||| i 2 (n) ||V(p||| z<4( n)X 

x (\\A<p%\\ LHn) + 2||An||U4 (a) ) (||n||Uoo (n) + M\\ L ~ { n)) 
+ l|VAv?fe|| L 2 (n) ||V^|| L 6 (n) x 

x (||Vpf |U 6(n) + 2||Vn||| i6(a) ) (||Vn||| L 6 (n) + ||V^||| L 6 (n) ) 

~ ll VA ^ll i2 (")II^H ff2 ( f2 )(ll n fcll« 2 ( f2 ) + II^IHff 2 (f2)) X 

x (||V^||U2 (n) + ||VA^||i2 (n) + HVny^n) + II VAn||| L2(n) ) 
+ II VA^||| L 2 (n) ||v|||ij2 (a) (||n||| H 2(a) + |bfe|| ff 2 (n) )x 

x (||A^||i2 ( n) + ||VA^|| L2(0) + ||An||| L 2 (a) + \\VAn%\\ L ^ a )) 
+ l|VA v ||| L 2 (n) ||^||| H 2 (a) (||n||| H 2(a) + ||v||| H 2 (n) ) 

< ||VA^|| La(n) ||^|| H 2 (n) (||n||| ff 2 (n) + ||^|| H2(n) )x 

x (MWhhq) + ||VA^||i2 (n) + \\n%\\ H 2 {Q) + ||VAn||| ia( n)) 

< (v + MWhhu) (l|nI||H2 ( n) + ||^||^(n))) l|VA^||| 2(n) 

H 2 (n) + ll«lll_H- 3 (o) + l|VAn||| L 2 (f2 )^ x 
x (\\n%\\ H Hn) + MWh^q)) , 
for all T] > 0, for some C, q > 0. 
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Also, for all rj > 0, there is C v > such that 

7J 2 , a = -(VA^I | 2Vn|AnS^) - (vA<p£ | |Vn|| 2 VpS) 

< l|VA ¥ j|||i a (n)(2||Vn|||i-(n)||An|||£a(n)||¥Jfc||z,-.(n) 

(420) +l|Vnf||£ 6(n) ||V^|| L6(n) ) 

< 77||VA^|||| 2(n) 

+ ^(||n|||^ 2(n) + ||VAn|||| 2(fi) + l)||n|||| 2(n) ||^|||| f2(n) . 



Summing up (4.19) and (4.20), we get: there is C > 0, and for all -q > 0, there 
is C v > (with C and C v depending on ||fto|U o °((o,oo),ii' 2 (n))) such that 

(4.21) 

^3<^(*7 + ||VfclU»(n)(l + ||vIllfl»(Q)))l|VA^||i, (n) 
Estimating Now, 

II 3 = 1/3,1+7/3,2, 

with 

// 3 ,i = -a(AV| I m% A (m| A M™l)) - n% A (n| A /i d (n|))) 
// 3>2 = -a(AVl I mf A (m| A /i ext (i)) - n% A (nf A WO))) • 

Concerning // 3j i, first write m| = n e k -\-ip e k , then integrate once by parts, so that 
J/31 takes the form of a 7 2 scalar product between VA<p| and a sum of terms 
V(afrc), where a, 6, c may be n| (or h^{n k )) or (or hd(<p k )), and at least one 
of them is n| (or h&{n e k )). Estimating each of a, 6, c and their gradients in 7 6 , 
one gets: for all r\ > 0, there is C, ; > such that 

(4.22) // 3 ,i < ?7||VA^||| (l + IKII?P(n) + IKII M\\ 

Then, split 7/3,2, 

//3,2 = (A VI I ™l A (m% A (/lextW - /lext(0))" 



L 2 (fi) 



- a (A V| I m| A (m| A /i cxt (0)) - n% A (n£ A h cxt (0) 



L 2 (Q) 



The second term is estimated as 7/3,1. The first one is split into a sum 
involving n| A (n| A (h cxt (t) — ft. cx t(0))^ , and products of h cxt (t) with two terms, 
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one of them being (pf., and the other, ipf, or n £ k . This leads to: for all 77 > 0, 
there is C v > (also depending on ft. ox t) such that 

4 ^ Ih,2 <*7l|VA^.||£ 2(n) + C„((l + Kllfp,^) + M\\ 2 H 2 (n) ) \Hf„ H n) 

+ *ll"'llllr=(n) lbIIU 2 (n) 



Finally, summing up (4.22) and (4.23), we have: for all 77 > 0, there is C v > 
(depending on \\no\\ L ^((o,oo),H 2 (Q))) sucn tnat 

J/ 3 <»j||VA^||| a(n) 

(4.24) (t ,2 s 
+ Cr,[[l + M\\ 2 H ^n)) yUHHV+tMUHn))- 

Estimating II4. Integrating once by parts, we get 

Ih = -(VA^I I V[P k , H0,-)}(n e )) L2m - 

Thus, for all 77 > 0, there exists C n > such that 

(4.25) ^4<»7l|VA^||i 3(n) +C IJ r| il , 
with 

r% 1 = ||V[P fe , J"(0, ■)}(n e )\\l 3(n) — -> in L°°(0,T) for all T > 0, with E fixed, 

v ; fc— >oo 

thanks to Lemma 4.1. 
4.1.4 Conclusion 

^From (4.18), (4.21), (4.24) and (4.25), we deduce that there is a constant 
C > (depending on ||^o||z°°((o,oo),.ff 3 (n)))> an d for all 77 > 0, there is C v > 
(depending on 77, ||rao||L<~((o,<x>),if3(n)), ||'iext|li«. w -i,oo and ll^eoctlli^wl- 00 )' 
such that 

~ (HArtlli,^,) 
2g +(« - C(r/ + lbllk 2 (n)(l + M\\m { n)))) ||VA^|||| 2(n) 

«S ^(||^||^ ( n)(l + ||</>IHff2(fi))(l + IKIlH» ( n) + HVAn|||| 2(n) ) 

+ %IIIh=(Q) +r| >:L ). 

Sum up (4.10) and (4.26), to get: there is a constant C > (depending on 
II«o||l~((o,oo),// 3 (o))), and for all 77 > 0, there is C v > (depending on 77, 
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II«o||l-((o : oo),h 2 (o)), l|/iext||£oo W i.°° and 119*^11^=0^'°°)' such that 
(4.27) 

" " l Vllll^n)) + (a-G(?7+||^|| ff 2 (0) (l + ||^||| H2(0) )))||VA^|||2 f2( ^ ) 



2dt 



^(11^111^(^(1 + MflapijXl + Ikllllp(Q) + ||VAn|||| 2(n) ) 

withff = 11^,^(0, -)](n e )||2 — ^ in L°°(0,T) for all T > 0, with e fixed. 

Now, apply the following Gronwall lemma (the proof of which is postponed 
to Section 6.3). 

Lemma 4.2. There is a constant K > (depending on uq and h ex t) such that, 
for all c £ (0,1/K), setting t s = celn(l/e), there is £0 = eo(a,c, K) such that 
(4.27) implies: 



Vee(0,e o ], 3k(e)eN\ Vfc ^ k(e) 

sup M\\ 2 HHn) < + K\\f%\\ L ^ o) e-^ e K ^ VAP ^ao,^. 



4.1.5 Passing to the limit k — > 00 

For each e € (0,£o] fixed, by Lemma 4.2, the sequence (ip k )k^w is bounded in 
L°°((0, t e ), H 2 (£l)). Equation (4.5) then implies that the sequence (d t <p%)k£n* 
is bounded in L°°((0, t e ), L 2 (Qj). Furthermore, (4.27) shows that (<p k )kef>i* is 
also bounded in L 2 ((0, t E ), H 3 (n)). Aubin's Lemma (see [2], [10]) then implies 
that there is a subsequence of (<p%)keN* converging in £ 2 ((0, t e ),H 2 (Q)) towards 
some (p £ . 

Up to a subsequence, we may assume that (dttp k )k£N* also converges weakly 
in L 2 ((0, < e ), L 2 {Q)) towards d t f £ ■ As k goes to 00, Pkn £ converges towards n £ 
inC([0,i £ ],iJ 2 (fi))ni/ 1 ((0,t e ),L 2 (r2)). Thus, (m|) fceN * converges towards some 
rrf in L 2 ((0, t E ), H 2 (il)), with (9 t m|) fceN » converging weakly in L 2 ((0, t s ), L 2 (il)) 
towards 9 t m e . This is enough to pass to the limit in (4.4), so that m £ is solution 
to (4.2). With e fixed, showing that m £ is continuous in time with values in H 2 
is standard, as well as uniqueness and stability properties: see [4], or [1]. 

Finally, passing to the limit in Lemma 4.2 yields: 

snp||^||^ ( a)<^e K||VA " o|l ^«o-)xn, ; 
[o,*.] K 

which we write 

(4.28) sup 11^11^(0) < K'e x - cK . 

[0,*e] 
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4.2 Second step: following the slow dynamics after t e 

iFiom the local-in-time existence result, we know that, for each e £ (0,Sq), 
there is t £ > t e such that m £ exists, as a solution to (4.2), in C([0, t e ], H 2 (Cl)) n 
i 2 ((0, t £ ), H 3 (Q)). We shall show, via a priori estimates, that t £ ^ T (possibly 
reducing eq). 

^From (3.1), and with T from (4.3), we deduce that, on [0,T] x f2, 

ed t m cq - aAm cq = F(t, m cq ) + ed t m cq . 

Substracting to (4.2), we get (on [0,t £ ] x fl): 
(4.29) 

( (ed t - aA)(m £ - m cq ) = (c(m eq ) + K(m cq )j (m £ - m oq ) + ed t m cqi 
I ^(m 5 - m eq )| sn = 0, 

and we consider the associated Cauchy problem with data given at time t s . The 
data at time t e = celn(l/e) satisfy (using (4.28) and (4.1)): 

(4.30) 

\\(m £ - m cq )(t e )\\ H 2 (n) sC \\(m £ - n £ )(t e )\\ H 2 {n) + \\n £ (t £ ) - m cq (0)\\ H 2 (n) 

+ ll^cq(O) - meq(t e )||ff»(n) 

s? K'e 1 -^ + ||no(cln(l/e)) - m eq (0)|| H2(n) 

+ ||Weq(0) - m eq (t e )||ij=(n) — >0. 

Here, for all S € H 2 (fl) and t G [0,T], 

£(t, m cq (t)) S =a|Vm eq (t)\ 2 5 + 2a(Vm cq (t) • V<j)m oq (t) 

+ <5 A /jr(t, meq(t)) + m eq (t) A (A5 + /i d (5) 
(4.31) -a<5A (m cq {t) A (/i d (m oq (t)) + /i cx t(*))) 

- am cq {t) A f<5 A f/l d (TO eq (i)) + ^ext(*) 

- am cq (() A (m cq (i) A /i d (5) 



and 
(4.32) 

n(t,m e q(t)) (6) = 2afVm eq (t) ■ VS^jd + a\VS\ 2 5 + S A (A<5 + /i d (5)) 

- a<5 A f<5 A (/i d (m eq (i)) + /lext(t) jj - aS A fm eq (t) A /i d (<5)J 

- am eq (i) A (s A /i d ( 5 )) - a5 A ( S A h d (S)j . 
In the sequel, we consider 

S £ := m £ - m eq G C([0, i £ ], ff 2 (ft)) n L 2 ((0, t £ ), # 3 (tt)), 
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and we simply prove that in (H 2 ) energy estimates, the term due to the resid- 
ual lZ(m cq )5 e is dominated by the terms due to aAS £ and to the linear term 
C(m eq )d e . We thus come back to the Galerkine approximation Sf, of 5 e , as 
in Paragraph 4.1.1. Take the L 2 (Vt) scalar product of the equations with 
(5| and A 2 Sf. and integrate by parts. Estimating (<5| | 7£(m eq )(<5|)) L 2 (n) is 
straightforward. Due to the continuity properties of h d on Sobolev spaces, 
(A 2 (5| | 72.(m oq )((5|)) i 2( ) produces three kinds of terms. Dropping the expo- 
nent e and subscript k (and using the notation L(v\, . . . , v n ) for any n-linear 
application), we examine each of them. 

From 5 A A h d (6)j . We have 

(A 2 <5 | L(S,S,S)) LHn) = (AS | AL(5,S,S)) LHn) s; \\A5\\ LHn) \\AL(S,S,5)\\ LHQ) , 

which is bounded from above by CH^Hlp^), since H 2 (Vt) is an algebra. 

In the same way, the terms of the form (A 2 <5 | L(S, (>))l 2 (q) are controlled by 
ll^ll^^) • This rules out the terms from SAh d (S), SA^SA(jid(m cq (t))+h cxt (t)j^j, 
S A (m cq (i) A h d (S)^j and m cq (t) A (s A h d (5)^j . 

From \VS\ 2 S. Write 
(A 2 6 | L(VS,VS,S)) L 2 (n) = 

- (VAS | L(VS, V6,V6)) LHQ ) - (VA<5 | L(A5, VS, 6)) L > [a) . 

Then, 

\(VA6 | L(V5,VS,VS)) L 2 {n} \ < \\VAS\\ LHQ) \\L(VS, VS, VS)\\ LHU) 

<C||VA*|| £3(tl) ||W||!a (n)1 

and by Sobolev's inequalities, ||V(5|| £ 6(q) is controlled by ||<5||ij2(m. 
Also, 

|(VA<5 | L(AS.yS,S)) L 2 (n) \ <: C\\VA5\\ L 2 {n) \\A6\\ L 2 in) \\V5\\ L ~ (n) \\5\\ L ~ {n) . 
^From the estimate 

l|V*|U- (n) < \\VS\\ L2{n) + \\VAS\\ L , m , 

we get 

\(VA8 | L(AS, V6,6)) L , m \ < C(||VA«|| x , (n) ||*||| P i (n) + ||VA5||| 3(n) ||*||^ (n) ). 
This leads to 

(A 2 6 | L(V6,V6,8)) L2m sc C(||VA<J|| L 2 (n) ||cJ|| 

H 2 {n) + l!^A(5|j 2 2(f2)ll^llff2(fj)). 

In the same way, we have 

(A 2 S | (Vm eq (i) • V5)S) < C\\VAS\\ L 2 m \\S\\ 2 H2{n) . 
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The 6 A AS term. Again, 

(A 2 S | L(S,AS)) LHn) = — (VA<5 | L(V6, A5)) L2(n) - (VAS \ L(5,VA5)) L * m , 
and as above, we get 

(A 2 S | L(5,AS)) L * m ^ C(\\VA5\\ L 2 {n) \\S\\ 2 H2{n) + ||VA<5||| 2(n) |H| H2(n) ). 
Finally, there is C > 0, and for all r\ > 0, there is C, ; > such that 

(S £ k \TZ(m cq )(8l)) H 2(n) ^(j) + C\\Sl\\ H 2^ n) + C\\S £ k \\ 2 H2 ^ n ) )||VA^||i, (n) 

+ C?)(ll^fellff2( a ) + ||<5fc||ff2( ) + \\8 e k\\H 2 (n))- 

Let fc go to infinity, so that the above estimate applies to S e instead of <5|, up 
to the local existence time t e obtained via the convergence of the Galerkine 
scheme. Coming back to (4.29), still with 5 e — m s — m cq , we get, using (2.5): 
there is C > 0, and for all rj > 0, there is C v > such that 

(4.33) 

2 S (l^ll W») + (<*-*- c(i + W\\ HHn) w\\ HHn) ) ||V^|| 2 H2(0) 

<(c„(l + W\\ HH n)) 3 W\\ H 2 (n) - C hn ) W\\ 2 H2m + e 2 \\d t m cq \\ 2 H2(n) . 

As in the proof of Lemma 4.2, fix i] E (0, a), and consider the time t e ^ t £ up 
to which, in (4.33), the parenthesis in front of || V<5 e ||^- 2 (Q) (resp. ||^ E ||#2(f2)) 
remains positive (resp. less than — Cn n /2). We have, for t £ (t e ,t £ ): 



Gronwall's lemma then implies that 



sup ||<5 e ||ff2 (a) \\6 e {t e )\\ 2 H2(n) +2 e Tsup ||9im eq ||^ 2(n) , 

lU.fi] [0,T] 

so that, for e small enough, we get t e ^ T, and sup[ te T ] — >Q. This 
finishes the proof of Theorem 2.1. □ 



5 Proof of Proposition 2.2 and Corollary 2.3 

Proof of Proposition 2.2. For any T > and n E C([0, T],H 2 (Q)), it is 
equivalent for n to be solution to (2.8) or to 

(d t - aA)(n - m cq ) = ^£(0, m cq ) + K(0, m cq )j (n ~ m cq ), 
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with the same initial and boundary conditions. The operators £(0, m eq ) and 
TZ(0,m cq ) from (4.31) and (4.32) do not depend on time, now. Arguing as in 
Section 4.2, we get an estimate analogue to (4.33), 

(5.1) 

2^(l|n-m eq ||| 2(n) ) 

+ (a-r)-C(l + \\n - m eq ||H 2 (fi))l|rc - m e q||jj2 (n) J ||V(n - m eq )||^2( n) 
4:\C V (1 + ||n - m oq ||,H-2 (rl) ) 3 ||n - m eq \\ H 2^ - C*n n )||n - m cq \\ 2 H2 ^ n y 

Once i] £ (0, a/2) is chosen, take r] > such that, when ||m — m eq || h 2 ((i) ^ Voi 
the parentheses in front of || V(n — TO e q )|| j^q) and in front of \\n — "m-cqW 2 ^^ 
are positive and negative at t = 0, respectively. The bootstrap argument then 
shows that n e C([0, oo), H 2 (tt)), and that n(t) converges in H 2 (Vl,S 2 ), as t 
goes to oo, towards m cq (to)- 



(5.2) \\n(t) - m cq \\ H 2( n) < r] e 



-ci 



for some C £ (0,Cii n ) depending on 770. Coming back to (5.1), we see also that 
V(n-m cq ) € L 2 ((0, 00), il 2 (n)). □ 

Proof of Corollary 2.3. When m oq (0) is constant over fl, Proposition 2.2 
ensures there exists some 770 > such that for all m £ ifjy(£2, S* 2 ) satisfying 

|| mo - m cq \\ H 2( n) ^ r/ , 

Asumption (ii) in Theorem 2.1 holds true. Furthermore, estimations (5.1) and 
(5.2) show that the corresponding function n has norms in L°°((0, 00), H 2 (£l)) 
and in L 2 ((0, 00) x f2)) controlled in terms of m oq (0) and 770 only. Thus, Eq in 
the proof of Theorem 2.1 may also be chosen depending on m eq (0) and 770 only, 
uniformly with respect to tuq. □ 



6 Appendix 

6.1 About the dissipation property (2.5): proof of Lem- 
mas 3.5 and 3.6 

Let S £ C([0,T],H°°(Q)) be such that |m eq + <5| = 1 and d v S\ ga = d v A5\ gn = 0. 

Then, 

(6.1) 

C(mf )6 = (\Td)SAu±uA(AS+h d (5)) + a(d=FA)uA((5Au) - auA(uAh d (S)). 
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6.1.1 L 2 estimates 

Take the L 2 (il) scalar product of (6.1) with S. This yields 



(£(m± ) 5 | 8) mu) = ±JJ-{uAAS)± J f-(uA h d (S)) 



'l 2 (o) 

(6.2) 

+ a(dTA) / \5Au\ 2 ~a I (u ■ h d (S)){u ■ S) + a I S ■ h d (S). 
Jn Jn Jn 

First consider the case of m£ q . Denoting n the exterior normal vector to f2, the 
first term in the right-hand side of (6.2) is equal to 



3 n 3 n n 

(6.3) V / S-di(uAdi6) = V / 6 ■ (u A di8)rn = / 5 ■ (u A d n 5) = 0. 
„■_ -, Jon Jon 



3 „ 3 

/ S-di(uAdiS) = Y 
i=1 Jn i=1 Jon J do, 

Since h d is continuous on L 2 with norm 1, the second term is bounded from 
above by H^H^^)- Similarly, due to the non-positivity of h d , the last term is 
non-positive. In the two other terms, we inject the identities 

(6.4) \5\ 2 = -2u-5 and \5 A u\ 2 = \5\ 2 - , 

which stem from the equality \u + S\ = 1. This leads to 

(C(mt q )S\S) L2in) ^\\S\\ 2 L2(n) +a(d-\) [ (\S\ 2 |<5| 4 /4) 

(6 ' 5) +f f(u-h d (S))\5\ 2 

= (l + a(d-A))||<y||i a(n) +0(||*||i a(n) ). 



In the case of m cq , we obtain in the same way 

(6.6) (£(m- ) 5 | 6) L2(n) ^ (a(X + d) - c) ||«5||| 2(n) + 0{\\5\\l* m l 
for some constant c depending on and a only. 

6.1.2 H 2 estimates 

Take the L 2 (f2) scalar product of the Laplacian of each term in (6.1) with A<5. 
This yields 

(A£(m±)6 | AS) L2{n) =(\Td) J AS ■ A{5 A u) ± J AS-A(uAAS) 

(6.7) ± f 5- A(uAh d {5)) + a{dTX) [ \AS/\u\ 2 

Jn Jn 



a AS ■ A(u A (u A h d (S))). 



24 



Since A(SAu) — (AS)Au = 0, the first term on the right-hand side vanishes. So 
does the second one, by the same argument as in (6.3). The equality |u + <5| = 1 
implies 

\AS A u\ 2 = \A5\ 2 - (|V<5| 2 - (6 ■ AS) 2 ) 2 , 

so that (6.7) gives, for : 
(6.8) 

(A£(m+) S | AS) L2(n) < -a(X - d)\\A6\\% m + c\\6\\ 2 H2(n) + O (||<5|| 3 ff2(0) ) , 

for some constant c depending on Q, and a only. Together with (6.5), we get 
finally 

(6.9) (£(m+ ) S | S) H2{n) sc - (a(X - d) - c) \\S\\ 2 H . 2{n) + O (||<5|| 3 ff2(0) ) , 



which concludes the proof of Lemma 3.5. □ 
In 

(6.10) 



In the case of m eq , we have 



(A£(m+)J|A£) i2(n) =- / A*.A(t»Ahd(*)) + a(A + d)l|A*||i 3(n) 

-a f AS ■ A(u A (u A h d (5))) + O (\\Sf H2{Q) 

•/fa 

which, together with (6.6), leads to Lemma 3.6. □ 



6.2 Proof of the commutator lemma 4.1 

Writting 

[P k , J-(0, -)](n) = (P fc - l)7-(0, n) + J"(0, n) - -F(0, P fc n), 

the result follows from the convergence of Pk towards 1 pointwise as an operator 
on i? 1 (n) (which rules out the term {P k — 1).F(0, n)) as well as on H^(£l), 
combined (to deal with .F(0, n) — -F(0, Ptn)) with the continuity of .F(0, •) from 
C([0, T],H 2 (fl)) n i 2 ((0, T), if 3 ) to L 2 ((0, T), if 1 ). 

The latter is a consequence of the continuity properties of and of Sobolev's 
cmbcddings, implying that H 2 (fl) is an algebra (so that all applications n i— > 
n A ha{n), n n- n A (n A /id(j'i)), i-> n A /i e xt(0), n n A (n A /i ex t(0)) 
are continuous on L°°((0, T), H 2 (fl) j), and that the product operation maps 
H 2 x H 1 to if 1 , so that n n- n A An and n i— > |Vn| 2 n are continuous from 
L°°((0,T),H 2 (n))nL 2 {{0,T),H 3 ) to L 2 ((0, T), if 1 ). □ 



6.3 Proof of Gronwall's lemma 4.2 

First, consider tef and e > fixed. Set <p E (t) = \\<p%\\ 2 H 2rn)' r ( t ) = 
iV (t) = ||VAP fc n (t)|| 2 2(n) and iV s (t) = ||VAn|(t)|| 2 2(0) , so that 

N £ (t) = N (t/e) and N a e i^Coo). 
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With C from (4.27), choose r\ G (0, a/{2C)). Hence, there exists k v £ (0,1) 
such that 

Vcp£[0,K v ], C(ri + ip(l + (p)) <a/2. 
Set K = 8C n (also from (4.27)), c £ (0, l/K) and t E = celn(l/e). Then, with 

t% = sup{£ £ [0,t e ] I Vi' G [0,t], </> £ (t') < k„} 

(i| > since </> e (0) = 0), we have: 

Vt £ [0, ecj) e '{t) < X((l + A £ (t))</> £ (t) + t + r(t)). 

^From this, we deduce: 

Vte[0,t£], ^ e (*)<^ ji 1 ' + r ^) exp {§ I ( 1 + N e (t"))dt"^\dt / 

(J* ^(t> + r ( t >fj cxp (^^) e K|Ar ° llll<0 -» 

(^ + flk|U 1( 0,T )) e ^"^"-«o.-, 



with T = ceohi(l/eo) ( an d £0 is chosen below). Now, since c € (0, l/K) and 
i e = celn(l/e), 

/ 1-cK \ 

V< G [0,min(t|,i £ )], e (i) < f i-^- +if||r|| L1(0 , To)£ - 1 -^J e *ll"oll^co,«o , 
which is less or equal to as soon as e belongs to (0, eo] j for 

/1 N l/(l-cJf) 

and k greater than K(e) such that 

Vfc > Jf( e ), ||f||U 1(0 ,T ) < ^K n £ -^e- K ^Ho^ 

(which is possible by Lemma 4.1). For this choice of e and fc, we thus have 
t| ^ t E , and the result follows. □ 
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